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We show how to calculate orbit implications (based on 
Nielsen- Thurston theory) for horseshoe type maps arising in 
chaotic time series data. This analysis is applied to data from 
the Belousov-Zhabotinskii reaction and allows us to (i) predict 
the existence of orbits of arbitrarily high periods from a finite 
amount of time series data, (ii) calculate a lower bound to the 
topological entropy, and (iii) establish a "topological model" 
of a system directly from an experimental time series. 

05.45.+b, 47.52.+j 



I. INTRODUCTION 

Braids arise as periodic orbits in dynamical systems 
modeled by three-dimensional flows jl]-{| . The existence 
of a single periodic orbit of a dynamical system can imply 
the coexistence of many other periodic orbits. The most 
well-known example of this phenomenon occurs in the 
field of one-dimensional dynamics and is described by 
Sarkovskii's Theorem ||. Less well-known is the fact 
that analogous results hold for two-dimensional systems 
||. In one-dimensional dynamics it is useful to study 
the period (or the permutation) of an orbit In two- 
dimensional systems it is useful to study the braid type of 
an orbit Q . Given this specification, we can ask whether 
or not the existence of a given braid (periodic orbit) forces 
the existence of another; as in the one-dimensional case, 
algorithms have recently been developed for answering 
this question 

As originally observed by Auerbach et. al., unstable 
periodic orbits are available in abundance from a single 
chaotic time series using the method of close recurrence 
[jlO|,[ll]]. By a "braid analysis" we propose to analyze 
a chaotic time series by first extracting an (incomplete) 
spectrum of periodic orbits, and second ordering the ex- 
tracted orbits according to their orbit forcing relation- 
ship. As shown in this paper, it is often possible to find 
a single periodic orbit (or a small collection of orbits) 
which forces many orbits in the observed spectrum. This 
single orbit also forces additional orbits of arbitrarily high 
period. This analysis is restricted to "low-dimensional" 
flows (roughly, flows which can be modeled by systems 
with one unstable Lyapunov exponent), however it has a 
strong predictive capability. 

We would also like to point out that this analysis gives 
us an effective and mathematically well defined "pruning 
procedure" for chaotic two-dimensional diffcomorphisms 



|12| . Instead of asking for rules describing which orbits 
are missing (pruned), we instead look for those orbits 
which must be present. All orbits of a low enough period 
(say up to period 10) appear to be succesfully perdicted 
by this method. This procedure will usually miss orbits 
of higher period, however from an experimental view- 
point the low period orbits are the most important and 
accessible. Orbits of low period often force an infinity of 
other orbits. This is illustrated in onc-dimcnsional dy- 
namics by the famous statement "period three implies 
chaos" [Of . An analogous statement in two-dimensional 
dynamics is that a non-well-ordered period three braid 
implies chaos [Q. 

This paper is organized as follows. In section II we 
briefly review the theory of orbit forcing in one- and two- 
dimensional maps (three-dimensional flows). In section 
III we show how the braid analysis method works by ap- 
plying it to times series data generated from the Rossler 
equations. This section also discusses some useful nu- 
merical and symbolic refinements to the method of close 
recurrence. In section IV we apply the braid analysis 
method to data from the Belousov-Zhabotinskii reaction. 
Our analysis builds directly on the original topological 
analysis of this data set due to Mindlin et. al. p5[ . In 
section V we offer some concluding remarks by indicating 
how this method can be extended to low-dissipation and 
conservative systems. In the examples studied in this pa- 
per we do have good control of the symbolics. In princi- 
ple, though, this method does not require good control of 
the symbolics (a good partition) and can thus overcome 
some of the current difficulties associated with finding 
good symbolic descriptions for (nonhyperbolic) strange 
attractors Jh|. 



II. ORBIT FORCING THEORY 

In this paper we study in detail the orbit structure of 
horseshoe type attractors. In particular, we will make 
a good deal of use of those features of one-dimensional 
dynamics which must always carry over to two dimen- 
sions. In our review we closely follow the point of view 
adopted in the Ph.D thesis of T. D. Hall Q. Thus, none 
of the results in this theoretical section are new. How- 
ever they are presented in such a way as to make them 
more useful for those wishing to apply these theoretical 
results to experimental data sets, as is done in the later 
sections of this paper. Let us begin by summarizing some 
well-established |3 18 definitions and results pertaining 
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to the periodic orbit structure of maps of the line. Let 
F : R — » R be a continuous map, and R = {x\, . . . , x n } 
be a period n orbit of F with its points labeled in such 
a way that x% < x-i < ■ ■ ■ < x n . Then there is a cyclic 
permutation 7r = tt(R, F) associated to R, defined by 
F(xi) =x 7r (i)- A relation < i on the set C of all cyclic per- 
mutations can be constructed as follows: if a, r G C, then 
t <x cr if and only if every map F : R — > R which has 
a periodic orbit R with tt(R, F) = a also has a periodic 
orbit 5 with tt(S, F) = r. If this is the case then we shall 
say that the permutation a forces the permutation r. The 
relation <i is a partial order which can be calculated, in 
the sense that there is a simple algorithm to determine 
whether or not r <j a for any two cyclic permutations r 
and cr: unfortunately this algorithm takes a long time to 
execute for orbits of high period, and it is difficult to use 
it to analyze the global structure of the partially ordered 
set (C, <i) 0. There is, however, an interesting subset 
of C onto which the restriction of <i is well- understood. 
Let UM C C be the set of unimodal permutations: that 
is, the set of permutations of periodic orbits of the tent 
map (or one-dimensional horseshoe) defined well into the 
parameter regime where a complete horeshoe exists (a 
full shift on two symbols) without prunning 

( 3x < x < 1/2 

F(x) = I 3(1 -x) 1/2 < x < 1 
I otherwise 

(equivalently, UM is the set of permutations which can be 
realized by periodic orbits of unimodal maps whose turn- 
ing point is a maximum). The partial order <i restricts 
to a linear order on UM which is described by kneading 
theory jllj. Indeed, this observation leads to a state- 
ment of topological universality in the one-dimensional 
context: there is a universal order in which the periodic 
orbits are built up in families of unimodal maps of the 
line, and this order is calculable using kneading theory. 

A periodic orbit R of a homeomorphism of the plane 
can be described by its braid type ; the reader who lacks 
the mathematical background necessary to appreciate 
the following definitions should simply regard the braid 
type as a two-dimensional analog of the permutation, and 
the term 'pseudo-Anosov' as one which, while essential 
for the accurate statement of theorem 1, is not essen- 
tial for applying the results found in later statements. 
Very roughly, braids can be divided into three types, the 
so-called finite order braids — braids which only force 
the existence of a finite number of other braids (periodic 
orbits), — pseudo-Anosov braids — braids which must 
force the existence of an infinite number of other braids, 
and reducible braids — braids which can be decomposed 
into distinct components which fall into one of the first 
two types. Pseudo-Anosov braids play an important role 
in the study of chaotic three-dimensional flows since, at 
the topological level, they force the coexistence of an in- 
finite number of other periodic orbits. The term pseudo- 
Anosov arises within the context of the Nielsen- Thurston 



classification of isotopy classes of surface automorphisms, 
which is the main tool used in the proofs of the results 
presented: the reader interested in this powerful theory 
could consult Ref . |2Q] , in which there is also an extensive 
bibliography. An exposition of the relevance of Nielsen- 
Thurston theory to two-dimensional dynamics, with par- 
ticular reference to braid types, can be found in Ref. [Q. 

Let / and g be orientation-preserving homeomor- 
phisms of the plane which have periodic orbits R and S 
respectively. We say that (i?, /) and (S, g) have the same 
braid type if there is an orientation-preserving homeomor- 
phism h : R 2 \ R — > R 2 \ S such that h o f\n2\ R o h^ 1 
is isotopic to g|R2\5. Having the same braid type is an 
equivalence relation on the set of pairs (i?, /), and we 
write bt(i?, /) for the equivalence class containing (R, /), 
the braid type of the periodic orbit R of /. We say that 
the braid type bt(i?, /) is of pseudo-Anosov, reducible, or 
finite order type, according to the Nielsen-Thurston type 
of the isotopy class of / : {S 2 ,R\J{oo}) -» (S 2 ,R\J{oo}) 
(see Ref. |2l[] for a statement of Thurston's classification 
theorem for homeomorphisms relative to a finite invari- 
ant set). By analogy with the one-dimensional theory, 
we define a relation <2 on the set BT of all braid types 
as follows: if (3, 7 S BT, then 7 <2 P if and only if every 
homeomorphism of the plane which has a periodic orbit 
of braid type (3 also has a periodic orbit of braid type 
7. If this is the case, we shall say that the braid type 
(3 forces the braid type 7. The relation < 2 is a partial 
order J2^|: as in the one-dimensional case, there is an 
algorithm for calculating this order (the so-called "train 
track algorithm") Q, but it is difficult to use it to ana- 
lyze the global structure of the partially ordered set (BT, 
<2). Still, using a version of this algorithm due to Bestv- 
ina and Handel 0, it is possible to calculate the orbits 
forced by an individual braid up to say period 10 or so 
by hand, and we have done this for several low period 
pseudo-Anosov horseshoes braids in order to find their 
associated topological entropies (see Table |). Also, an 
automated version of this algorithm promises to extend 
these calculations to higher period orbits p3[ . 

We restrict our attention to the subset HS C BT of 
horseshoe braid types: that is, those which are realized 
by periodic orbits of the horseshoe map / |24|] . The de- 
termination of constraints on the order in which periodic 
orbits can be built up in the creation of horseshoes is 
tantamount to the analysis of the partially ordered set 
(HS, <2); and an important observation in discussing 
the topological universality of two-dimensional systems 
is that, unlike (UM, <i), the partially ordered set (HS, 
<2) is not linearly ordered. 

The connection between the ordered sets (UM, <i) 
and (HS, <2) is that periodic orbits of both the tent 
map and the horseshoe have a symbolic description. We 
associate to points in the non- wandering set of each map 
an infinite sequence of 0's and l's in the usual manner, 
and describe a period n orbit R by its code cr, given 
by the first n symbols of the sequence associated to the 
rightmost point of R. We shall identify periodic orbits of 
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the two maps which have the same code, using the same 
symbol R to denote each orbit. Such an identified pair R 
will be referred to as a horseshoe orbit (or simply as an 
orbit): the fact that the orbit is periodic will always be 
assumed. We shall write cf? for the semi-infinite sequence 
crCrCr ■ ■ ■ (which is associated to the rightmost point of 
R, regarded as a periodic orbit of the tent map); and 
00 c^? for the bi-infmite sequence . . . crcrCr . . . (which is 
associated to the rightmost point R, regarded as a peri- 
odic orbit of the horseshoe). 

It will be convenient for us to make another identifica- 
tion. If R is a period n orbit, then denote by cr the code 
obtained by changing the last symbol of cr. If cr is not 
the two-fold repetition of a sequence of length n/2, then 
it is the code of a period n orbit R' which can be shown 
to have the same permutation and braid type as R. We 
identify R and R 1 in what follows, using the same sym- 
bol R to denote both, referring to them as a single orbit. 
The code cr of the identified pair will always be chosen to 
end with a zero (this is to ensure that the sequence 
contains the group of symbols . . . 010 . . ., which is nec- 
essary for the algorithms which follow). After making 
this identification, distinct horseshoe orbits have distinct 
permutations: however, it is possible for several different 
orbits to have the same braid type. With this in mind, 
we say that a function defined on the set of all horseshoe 
orbits is a braid type invariant if it takes the same value 
on any two orbits of the same braid type. 

Given two horseshoe orbits R and S, we shall write 
S <i R if and only if n(S, F) <t ir(R, F); and we shall 
write S < 2 R if and only if bt(5, /) < 2 ht(R, /). We re- 
fer to <i and <2 as the one- and two-dimensional forcing 
orders respectively. Our aim is to analyze <2 using the 
well-understood properties of <i. Notice that the orders 
<i and < 2 cannot strictly be regarded as being defined 
on the same set: if R and S are two horseshoe orbits 
which have the same braid type but different permuta- 
tions, then they must be regarded as being equal when 
considering the order <2, but as being distinct when con- 
sidering the orbit <i. For example, the orbits /3X2 and 
,Sg in Table | are two such orbits (for an example analo- 
gous example in the orientation-reversing Henon map see 
Ref. H). 

We say that a horseshoe orbit R is quasi-one- 
dimensional (qod) if for all orbits S we have R >± S ==>■ 
R > 2 S. As proved in Ref. [jl7|, the following result 
provides a simple classification of the qod orbits which 
have pseudo-Anosov braid type. Let Q denote the set 
Q n (0, 1/2) of rationals lying (strictly) between and 
1/2. Given such a rational q — m/n (in lowest terms), 
we write P q for the period n+2 orbit which has code c q = 
10 « 1 ( g ) 1 2 « 2 ( 9 ) 1 2 _ _ _ i2 "m(g) 10) where Ki(q) = [1/qj - 1, 
and Ki{q) = [i/q\ — [{i — l)/<zj — 2 for 2 < i < m (here 
denotes the greatest integer which does not exceed 
x). Thus, for example, c 2 /5 = 1011010. Then we have 
the following result of T. D. Hall (rj: 



Theorem 1 q < — ► P q is a one-to-one correspondence be- 
tween Q and the set of quasi-one-dimcnsional horseshoe 
orbits which have pseudo-Anosov braid type. Moreover 
P q >2 Pq' ^=>q' >q for all q, q' € Q. 

It is a consequence of quasi-one-dimensionality that 
P q >2 P q ' •<=>■ Pq >i P q ' , and therefore the braid types 
of the orbits P q form a linearly ordered subset of BT. The 
final statement of the theorem tells us that the ordering 
within this subset is simply the reverse of the usual or- 
dering on Q. 

Theorem 1 allows us to quickly identify a very useful 
subset of pseudo-Anosov braids and a program is eas- 
ily written to quickly generate all the horseshoe qod or- 
bits of any desired period. Moreover, Theorem 1 can be 
used to define an invariant of horseshoe braid type: the 
height q(R) of a horseshoe orbit R is the unique clement 
of [0, 1/2] with the property that for all q G Q we have 
q < q(R) P q >i R and q > q(R) => P q ^1 R (thus 
q(R) — sup{g G Q : P q >i R}). Because the orbits P q 
are qod, it follows immediately that q < q(R) => P q > 2 
R. In fact, it can be shown that the second property in 
the definition of height also holds for the two-dimensional 
forcing order: that is |i~7| ], 

Theorem 2 Let R be a horseshoe orbit and q G Q. 
Then q < q(R) => P q > 2 R and q > q(R) => P q ^ 2 R- 
In particular, height is a braid type invariant. 

We can also find the finite order braids because the finite 
order braids of the horseshoe are exactly those whose 
period equals the denominator of the height |]l7| , p6| . 
Because the height is defined in terms of the one- 
dimensional order, it can be calculated using kneading 
theory. We first define the height q(c) G (0, 1/2] of a semi- 
infinite sequence which begins c = 10 . . . , and which con- 
tains the group of symbols . . . 010 . . .. To do this, write 
c in the form c = 10 Kl l^i0 K2 1 M2 . . ., where m > and //, 
is either 1 or 2 for each i, with /ij = 1 only if ftj+i > 
(thus the Ki and \ii are uniquely determined by c, and 
the fact that c contains the group . . . 010 . . . means that 
\i s = 1 for some s). Let 

Ir[C) = Gr + EU^' (2r-l) + Er=i«i. 

for each r > 1, and let s be the least positive inte- 
ger such that either fi s = 1 or n^ x 7i(c) = 0. Then 
nf = iii(c) = (x, y] for some x and y: we define q(c) = y if 
fi s = 1 or J s _|_i(c) > n| =1 (c), and q(c) = x if jji s = 2 and 

7 s+ i(c) < n? =1 i 4 (c). 

Theorem 3 q(R) — g(cg): in particular, q(R) is posi- 
tive and rational [0 . 

For example, let R be the period 20 orbit with code 
cr = 10000110001100001100. We have ki = 4, fix = 2, 
«2 = 3, fj,2 = 2, K3 = 4, /X3 = 2, K4 = 2, and ^4 = 1. Thus 
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h = (1/6,1/5], I 2 = (2/11,2/10], h = (3/17,3/16], 
and h = (4/21,4/20]. Since 4/21 > 3/16 we have 
Ix H h n h n h = 0, and h > I x n 1% n I 3 . Therefore 
= max(/i n h D J 3 ) = 3/16. 
There is a relationship between height and a well- 
established braid type invariant: the height of a horse- 
shoe is equal to the left hand endpoint of its rotation 
interval. A proof of this is given in Ref. fI7|| , where a 
practical algorithm for determining the rotation inter- 
val of a horseshoe is described. It can also be shown 
using the algorithm for determining the height that for 
each q — m/n (in lowest terms), P q is the only period 
n + 2 orbit of height q; thus it is the only orbit of its 
braid type. This observation enables us to determine 
exactly which of the orbits P q are forced by a given or- 
bit R in the two-dimensional order: we now present an 
algorithm for this purpose. If R is a horseshoe orbit 
then we define the depth r(R) e (0,1/2] n Q of R as 
follows: consider all groups of the form ... OHIO .. . or 
. . . 01010 ... in the sequence °° eg 3 ; suppose that there are 
I such groups g\,...,gi contained in one period of 00 eg 3 . 
If I = then r(R) = 1/2. Otherwise, for each i < Z let /j 
be the code obtained by starting at the last 1 in gi and 
moving forwards through °° eg 3 , and hi be that obtained 
by starting at the first 1 and moving backwards. Then 
r(R) = mini<i<jmax(g(/ i ),g(6 i )). 

Theorem 4 r(R) is the unique element of [0, 1/2] with 
the property that for all q € Q we have r(R) < q ==>• 
R >2 P q and r(R) > q => R~^ 2 P q . In particular, depth 
is a braid type invariant |17j] . 

The following corollaries from theorems 2 and 4 will 
be very useful in our braid analysis since they allow 
to us to apply one-dimensional theory to calculate two- 
dimensional orbit forcings. 

Corollary 1 Let R and S be horseshoe orbits. If 
r(R) < q(S) then R > 2 S. On the other hand, if 
q(R) < q(S) and r(S) < r(R) then R^ 2 S and S^ 2 R: 
thus orbits of the braids types R and S can exist inde- 
pendently of each other. 

In fact, a much stronger result can be proved: if r(R) < 
q(S) then every homeomorphism of the plane which has 
a periodic orbit of braid type R has at least as many peri- 
odic orbits of the braid type S as does the horseshoe (the 
corollary only says that every such homeomorphism has 
at least one such orbit). Thus all of the periodic orbits 
of the braid type S must be created before any of the 
periodic orbits of braid type R in any family of home- 
omorphisms leading to the creation of a horseshoe. We 
can use the following corollary to locate some of these 
orbits. 

Corollary 2 If a homeomorphism of the disc has a qod 
orbit Pq, and R is another orbit which is forced by P q in 



the one-dimensional order (which means that its height 
is > q) , then the homeomorphism must have at least one 
orbit of the braid type of R. If in addition the height 
of R is strictly greater than q, then the homeomorphism 
must have at least as many orbits of the braid type of R 
as does the horseshoe p7| . 

For example, if a period 7 orbit with the braid type of 
1011010 (height 2/5) is extracted, then there must be 
at least one orbit of the braid type of 10110 (since this 
braid type has height 2/5), and at least as many of the 
braid type of 101110 as in the horseshoe (i.e., 2) need 
exist (since the height of these orbits is also 2/5). On 
the other hand, both of the orbits in the period 6 pair 
101115 must exist since the height of these orbits are 1/2. 
A listing of the low period qod orbits and some of the low 
period orbits which they force is found in Table [n[ 

Since the qod orbits essentially inherit the one- 
dimensional forcing order we can use one-dimensional 
methods to calculate their topological entropies. This 
information can be used to obtain lower bounds for 
the topological entropy of a partially-formed horseshoe. 
More precisely, given a horseshoe orbit R, let h(R) de- 
note the smallest possible entropy of a homeomorphism 
having an orbit of the same braid type as R. Now we 
know that for all q > r(R) we have R > 2 P q , and hence 
the topological entropy h(R) > h(P q ). However, since P q 
is a qod orbit, it can be shown that h(P q ) is equal to the 
entropy of P q regarded as an orbit in one dimension: this 
can be calculated using standard transition matrix tech- 
niques. A Mathematica program based on the method of 
Block et. al. ]^8| is available from T. D. Hall which allows 
us to calculate the topological entropies of the qod orbits. 
A graph of the function (q,h(P q )) appears in Ref. |JL7| ; 
it is monotonic and discontinuous everywhere. Thus, we 
can go directly from a qod orbit to a lower bound for the 
topological entropy. For example, the orbit R with code 
cr = 10011010 has r(R) = 1/3: any partially formed 
horseshoe which includes an orbit of the braid type of R 
has entropy greater than h(P 2 / 5 ) ~ 0.442. By taking ra- 
tionals closer and closer to 1/3 we get better and better 
estimates (e.g., 5/14 gives 0.481). We can also use the 
pseudo-Anosov orbits which are not qod to get estimates 
(e.g., the orbit sf mentioned above has h ~ 0.498), how- 
ever in this instance we must compute the entropy from 
a train track and this can be a difficult calculation. 

In Table | we collect together some useful facts about 
braids in the horseshoe up to period eight. All the re- 
ducible orbits up to period nine only have finite order 
components. This is because the lowest period with a 
pseudo-Anosov orbit is five, and numbers less than or 
equal to nine do not have proper factors greater than 
or equal to five. As mentioned previously, the topologi- 
cal entropy for the pseudo-Anosov orbits which are not 
quasi-one-dimensional is calculated by finding the "train 
track" using the method of Bestvina and Handel JtJ . This 
method results not only in a topological entropy, but also 
an explicit Markov partition (a 'topological model') and 
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associated Perron- Frobenius matrix which is useful for lo- 
cating in phase space where the predicted periodic orbits 
are to be found. 



III. ROSSLER BRAID ANALYSIS 

A braid analysis of a low dimensional chaotic time 
series consists of four steps once an appropriate three- 
dimensional space is created (i) the periodic orbits 
are extracted by the method of close recurrence [^9],[30| , 
(ii) the braid type of each periodic orbit is identified and 
the orbits are ordered by their two-dimensional forcing 
relationship (iii) a subset of braids are selected 

which have maximal forcing and which force the orbits 
extracted in step (i), and (iv) if possible, an attempt is 
made to verify that some of the predicted orbits (not 
originally extracted in step (i)) are found in the system. 

In practice, steps (i) and (ii) are greatly simplified if 
the template or knot-holder organizing the flow can be 
identified using the procedure described by Mindlin and 
co-workers ||,[ll||3l],[32| . Knowledge of the template helps 
in obtaining the symbolic names of the periodic orbits 
and in calculating the forcing relationship for the specific 
braids in that template. For instance, if the template is 
identified as a two-branch horseshoe knot holder (as are 
all the examples studied in this paper), then the theory 
of qod orbits of section II can be applied to simplify the 
analysis. 

Although template identification is very valuable, it is 
not essential for a braid analysis. Nor is the symbolic 
identification of the extracted orbits. In the worst case 
a braid analysis does require that the the braid conju- 
gacy class of each extracted periodic orbit is identified 
(see Elrifai and Morton p3[ , or Jaquemard |54j for algo- 
rithms), and that the minimal Markov model (a 'train 
track' in the language of Thurston) can be constructed 
for each braid (see Bestvina and Handel M , Los Q , and 
Franks and Misiurewicz H for algorithms). Algorithms 
exist for both of these steps, although the most computa- 
tionally efficient version of the braid conjugacy algorithm 
is probably not an effective solution beyond Bg. 

To illustrate the braid analysis we consider a chaotic 
attractor of the Rossler equations, 

x = -(y + z) 

y = x + ey 

z = / + xz — \iz 

with e = 0.17, / = 0.4 and /i = 0.85. The Rossler equa- 
tions are integrated through 10 5 cycles and the return 
map is examined at the half plane X = {(x,y, z) : x < 
& y = 0}. The template is easily identified as a horse- 
shoe with zero global torsion. This template identifica- 
tion is verified by calculating the relative rotation rates 
and linking numbers of the extracted periodic orbits as 
described by Mindlin et. al. pIMpl. 



To extract the (surrogate) periodic orbits by the 
method of close recurrence we first convert the return 
map from the sequence of values {x n ,z n ) directly into 
the symbol sequence of 0's and Vs. In this particular in- 
stance, since the map is close to one-dimensional, a good 
symbolic partition is obtained by examining the maxi- 
mum value of the next return map formed from the pro- 
jection on the x-coordinate — (— x n , — x n +i) — at the 
surface of section. Orbits passing to the left of the maxi- 
mum are labeled zero, and those to the right arc labeled 
one. Next we search this symbolic encoding for each and 
every periodic symbol string. Every time a periodic sym- 
bol string is found we calculate its recurrence and then 
save the instance of the orbit with the best recurrence. 
For instance, in searching for the period three orbit '100' 
we search the symbolic encoding of the return map for 
any instance of '100' and its cyclic permutations '010' 
and '001', and every time this symbol string is found, we 
next calculate its recurrence, which for this period three 
orbit is e y=0 = (x n+3 - x n ) 2 + {z n+3 - z n ) 2 , and then 
save the orbit with the minimum e. The advantage of 
this procedure of orbit extraction is that it is exhaustive. 
We search for every possible orbit up to a given period. 
In these studies we searched for all orbits between periods 
1 and 16. 

The resulting spectrum of periodic orbits up to period 
eight is shown in Table [II . The orbits which are present 
in (the full shift) complete hyperbolic system, and not 
present in Table III, are said to be pruned. Our goal is 
to predict as best as possible the pruned spectrum from 
the chaotic time series. 

Before we discuss the braid analysis, though, it is in- 
teresting to consider the number of orbits extracted as 
a function of the number of points in the return map. 
As expected, we find that the number of orbits that can 
be extracted increases with the number of points in the 
return map. More importantly, our numerical results 
strongly suggest that using the method of close recur- 
rence it is possible to obtain all the low period periodic 
orbits embedded within the strange attractor after exam- 
ining only a finite number of data points, in this instance 
10 5 . For instance, after 10 4 points are examined our nu- 
merical results show that no new periodic orbits are found 
below period six. Similarly, after examining 10 5 points 
we believe we have found all orbits up to period eight. 
These results also caution us when working with small 
data sets — the extracted orbit spectrum is expected to 
miss orbits either because the orbit is pruned (it is not 
in the strange set) or because the sample of the strange 
set we are examining fails in providing a close enough 
coverage over the whole attractor. 

Using the results in section II, the finite-order and 
quasi-one-dimensional orbits in the extracted spectrum 
are easily identified from their orbit codes (see Table 
III). Not unexpectedly, we see a sequence of quasi- 
one-dimensional orbits of increasing entropy (decreasing 
height) — the maximal qod orbit is the period 16 orbit 
1011011011011010 with entropy h « 0.480804. All orbits 
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forced by this period 16 orbit up to period 8 are present, 
and none are missing. So this period 16 qod orbit already 
gives us a very good hyperbolic set with which to approx- 
imate to our (possibly nonhyperbolic) chaotic attractor. 
Can we do better? Doing better in this instance means 
identifying a pseudo-Anosov orbit which is not qod, but 
perhaps implies the maximal qod orbit found. Indeed, 
in this data set there is such an orbit, it is the period 8 
orbit 10010100 with entropy h w 0.498093. Again, the 
spectrum of orbits forced by this maximal pA orbit are 
consistent with the extracted spectrum which was exam- 
ined up to period 16. 

This data set is close to one-dimensional so kneading 
theory also does quite well for predicting the low period 
orbits. For instance, we could consider the period 3 orbit 
100, and this orbit (based on one-dimensional unimodal 
theory) also accurately predicts most of the extracted 
spectrum. However, this period 3 orbit is finite order, 
and we know from the results of Holmes and Whitley that 
there will be many (possibly high period) orbits which 
are forced by <i but not by < 2 |35|] (in fact, we know 
that in 2-dimensions 100 forces only itself and a fixed 
point, and in 1-dimension it forces orbits of arbitrarily 
high period). Thus, although one-dimensional theory is 
a useful guide in this instance to the low period orbits, it 
can not be safely applied to make predictions about high 
period orbits. 



IV. BELOUSOV-ZHABOTINSKII REACTION 
BRAID ANALYSIS 

To further illustrate the braid analysis method we ob- 
tained data from the Belousov-Zhabotinskii chemical re- 
action |56) . This is the same data set analyzed by Mindlin 
et. al. and consists of 65,000 equally spaced points which 
measure the time dependence of the bromide ion con- 
centration in the stirred chemical reactor. Following the 
techniques described by Mindlin et. al. we also embed- 
ded the scalar time series x(i), i = 1, 2, . . . , N in R 3 via 
a differential phase space embedding described by 



yx(i) = x(i) 
y 2 (i) = x(i) 
y 3 (i) = x(i) 



A*yi(i-1), A = 0.995 
x(i - 1) 



from which we reproduced the 3-dimensional attractor 
and return map shown in Figs. 5 and 6(a) of Ref. |jT5| . 
The attractor is a zero global torsion horseshoe. There 
are approximately 125 points per cycle so the return map 
consists of about 520 points. 

Our technique for extracting (surrogate) periodic or- 
bits from this time series differs somewhat from that 
described in Ref. p5| . We use the same procedure de- 
scribed for the Rossler data: first the return map data 
is converted into a symbol sequence of 0's and l's de- 
pending on whether the orbit passes to the left or right 
of the maximum value of the return map, and second an 



exhaustive search is performed for all possible periodic 
orbits between periods 1 and 15. Again, this data set 
is almost one-dimensional and the simple symbolic pre- 
scription just described leads to a unique and consistent 
encoding of all the periodic orbits we are able to extract. 
Since the data set (at the return map) is small, we choose 
not to pick an arbitrary cut off for e, the close return cri- 
terion. Rather, we report the best e we are able to extract 
for a given periodic orbit (see Table IV). By including 
this additional piece of information we can make a more 
selective judgement about which close returns are, and 
are not, good surrogates for periodic orbits. In this way 
we are able to locate a few more periodic orbits than 
were originally reported in Ref. p5|] . Also, one perhaps 
surprising result comes out of this extraction method. 
It is not uncommon to find orbits of high period with 
very small recurrences. For example the period 13 orbit 
1011011101110 has a recurrence of e = 0.000712 which 
is significantly better than almost all of the lower period 
orbits. 

A list of the extracted orbits and their Thurston types 



is presented in Table IV . As expected, there is a sequence 



of quasi-one-dimensional orbits of increasing entropy, the 
largest of which is the period 16 orbit 1011011011011011 
with entropy h rs 0.48084. In this particular instance, 
the period 16 qod orbit is in fact the maximal pseudo- 
Anosov orbit in the data set and it forces all the extracted 
orbits in this data set except for the finite order period 
three orbit 101. Indeed, a careful analysis of this data set 
suggests that the signal is subject to a small parametric 
drift which carries it between the strange attractor and 
a stable period 3 orbit (whose 1-dimensional entropy is 
h w 0.4812) @. 

Table [v| shows the number of predicted and extracted 
orbits as a function of period. The number of forced 
orbits which are not in the extracted data set increases 
with the period. As with the Rossler data set, we be- 
lieve the forced orbits which are missing could actually 
be extracted from the data set if we were given a longer 
time series which provides a better coverage of the entire 
attractor. 



V. CONCLUSION 

In retrospect, we find it remarkable that such a small 
subset of periodic orbits (which are rather easy to get 
from experiments) contain so much topological and dy- 
namical information about a (low-dimensional) flow. As 
previously demonstrated, a few low period orbits are suf- 
ficient to determine the template describing the stretch- 
ing and folding of the strange set fl p^|j3l| . The template 
provides an upper bound to the topological entropy and 
is, in a sense, a maximally (i.e., a full shift) hyperbolic 
set which we formally associate to a (possibly nonhy- 
perbolic) strange set p9|| . In this paper we show how 
a sequence of periodic orbits (and their associated hy- 
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perbolic sets) can be used to obtain a collection of finer 
and finer approximations to a strange set which is prob- 
ably not hyperbolic. Formally we might say that the 
hyperbolic set associated with each pseudo-Anosov braid 
is embedded within the strange attractor we are trying 
to describe in the sense that the (possibly nonhyperbolic) 
strange set must contain at least all the orbits forced by 
the extracted pseudo-Anosov braid. We indirectly dis- 
cussed two measures of the goodness of this approxima- 
tion — the difference in the topological entropy and the 
difference between the forced and extracted low period 
orbits. Using cither measure, we have seen that it is pos- 
sible to select moderate period orbits (say < period 20) 
which provide good hyperbolic sets with which we can 
approximate a (possibly nonhyperbolic) strange set. 

The dynamical information derived from an orbit de- 
pends only on its braid type. As mentioned in section 
III, the braid type of an orbit can be determined without 
obtaining a good symbolic description of the orbits in the 
flow. We will illustrate these techniques in a future paper 
in which we consider a braid analysis of the bouncing ball 
system in a low-dissipation regime [ |38[ . Our strategy for 
handling the cases where a good symbolic description is 
not easy to obtain is to find a complete set of braid invari- 
ants on the small set of braids of interest. For instance, 
as mentioned in Table |. the exponent sum (simply the 
sum of the crossings) is a complete braid type invariant 
for horseshoe braids up to period 8. 

We would also like to point out that it is common for 
a collection of finite order braids to be pseudo-Anosov. 
This suggests an alternative strategy: instead finding a 
single orbit with maximal implications, it should also be 
possible to find a collection of orbits (possibly with very 
low period) which force all the observed orbits and pro- 
vides yet another hyperbolic approximating set. This is 
also the subject of future investigations. 
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TABLE I. Low period horseshoe orbits. The exponent sum is a complete braid type invariant 
for all horseshoe orbits up to period eight (see Ref. [17] for the explicit conjugacy). The braid 
name notation is from Ref. [35]. Thurston braid types: finite order (fo), reducible (red), and 
pseudo-Anosov (pA), and quasi-one-dimensional (qod). The "red, fo" orbits are reducible Thurston 
types with finite order components. 



r 


cp 


TYp 


Type 


es 


q{ p ) 


h 


braid name 


symbolic name 


permutation 


Thurston type 


exponent sum 


height 


topological entropy 


Sl 


1 


(1) 


fo 





1/2 





S2 


10 


(12) 


fo 


1 


1/2 





S3 


10? 


(123) 


fo 


2 


1/3 





f'2x'2 


1011 


(1324) 


red, fo 


5 


1/2 





S4 


100? 


(1234) 


fo 


3 


1/4 





si 

si 

SB 


1011? 
1001? 
1000? 


(13425) 
(12435) 
(12345) 


fo 

pA, qod 
fo 


8 
6 
4 


2/5 
1/3 
1/5 




0.543535 



4 


10111? 


(143526) 


red, fo 


13 


1/2 





/3x2 

4 

si 


100101 
10011? 
10001? 


(135246) 
(124536) 
(123546) 


red, fo 
red, fo 
pA, qod 


9 
9 
7 


1/3 
1/3 
1/4 





0.632974 


si 


10000? 


(123456) 


fo 


5 


1/6 





4 


101111? 


(1453627) 


fo 


18 


3/7 





4 


101101? 
100101? 


(1462537) 
(1362547) 


pA, qod 
P A 


16 
14 


2/5 
1/3 


0.442138 
0.476818 


4 


100111? 


(1254637) 


pA 


14 


1/3 


0.476818 


4 


100110? 


(1356247) 


fo 


12 


2/7 





4 


100010? 


(1246357) 


P A 


10 


1/4 


0.382245 


4 
4 


100011? 
100001? 


(1235647) 
(1234657) 


P A 
pA, qod 


10 

8 


1/4 
1/5 


0.382245 
0.666213 


4 


100000? 


(1234567) 


fo 


6 


1/7 





/22 X 2 

4 

sl 


10111010 
1011111? 
1011011? 


(15472638) 
(15463728) 
(14725638) 


red, fo 
red, fo 
fo 


23 
25 
21 


1/2 
1/2 

3/8 







4 


1001011? 


(13725648) 


pA 


19 


1/3 


0.346034 


4 
4 


1001010? 
1001110? 


(13647258) 
(13657248) 


pA 
pA 


17 
17 


1/3 
1/3 


0.498093 
0.498093 


„6 


i nm i i i 
lUUlllli 


(12564(58) 


pA 


19 


1/3 


a o a a no a 
0.34D034 


4 


1001101? 


(12573648) 


pA 


17 


1/3 


0.498093 


/2x4 


10001001 


(13572468) 


red, fo 


13 


1/4 





4 


1000101? 


(12473658) 


pA 


15 


1/4 


0.568666 


4 


1000111? 


(12365748) 


P A 


15 


1/4 


0.568666 


4° 


1000110? 


(12467358) 


red, fo 


13 


1/4 





4 1 

„12 

ss 

« 13 

s 8 

„14 
Sg 


1000010? 
1000011? 
1000001? 
1000000? 


(12357468) 
(12346758) 
(12345768) 
(12345678) 


pA 
P A 
pA, qod 
fo 


11 
11 

9 
7 


1/5 
1/5 
1/6 
1/8 


0.458911 
0.458911 
0.680477 
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TABLE II. A sequence of quasi-one-dimensional (qod) or- 
bits with increasing forcing implications (entropy). Thus, ev- 
ery orbited forced by the period 7 qod orbit is also forced by 
the period 10 orbit, and so on. The notation "(*)" after an 
orbit indicates that one (but perhaps not both) of the sad- 
dle-node pair is forced. Both saddle-node partners will be 
forced by the next orbit in the qod sequence. 



Period 


Forced by (7) 
101101? 


Forced by (10) 
101101101? 


Forced by (13) 
101101101101? 


1 








2 


10 






3 








4 


1011 















6 


10111? 






7 


101111? 


101101? 




8 


10111010 
1011111? 


1011011?(*) 




9 


10111111? 

mn nmO 

lOllllOli 
10110101? (*) 


10110111? 




i n 
1U 


iniiiiiii u 

lUlllllllx 

101111101? 
101110101? 
101101111? (*) 


ini imm i i 
1U1 1U1U1 1 1 

101101101? 


ini ini ini 1 ' 


11 


1011111111? 
1011111101? 
1011110111? 
1011110101? 
1011011111? (*) 
1011011101? (*) 


1011010111? 
1011010101? 


1011011011? (*) 
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TABLE III. Periodic orbits extracted from time series data 
of the Rossler equations. All extracted orbits up to pe- 
riod 8 are shown. Between periods 9 and 16 only the ex- 
tracted quasi-one-dimensional (qod) are shown. The maximal 
pseudo-Anosov orbit (IOOIOIOi) is not qod, but it forces all 
qod orbits with height greater than 1/3. Thurston Types: 
finite order (fo), reducible (red), and pseudo-Anosov (pA). 



Period 


Name 


Height 


Entropy 


Type 


1 


1 


1/2 





fo 


2 


10 


1/3 





fo 


3 


101 


1/3 





fo 


3 


100 


1/3 





fo 


4 


1011 


1/2 





red 


5 


10111 


2/5 





fo 


5 


10110 


2/5 





fo 


6 


101110 


1/2 





red 


6 


101111 


1/2 





red 


6 


100101 


1/3 





red 


7 


1011111 


3/7 





fo 


7 


1011110 


3/7 





fo 


7 


1011010 


2/5 


0.442138 


pA, qod 


7 


1011011 


2/5 


0.442138 


pA, qod 


7 


1001011 


1/3 


0.476818 


pA 


7 


1001010 


1/3 


0.476818 


pA 


8 


10111010 


1/2 





red 


8 


10111110 


1/2 





red 


8 


10111111 


1/2 





red 


8 


10110111 


3/8 





fo 


8 


10110110 


3/8 





fo 


8 


10010110 


1/3 


0.346034 


pA 


8 


10010111 


1/3 


0.346034 


pA 


g 


1 001 01 01 


1 /3 


498093 


nA 


8 


10010100 


1/3 


0.498093 


pA 


9 


101111010 


3/7 


0.397081 


pA, qod 


9 


101111011 


3/7 


0.397081 


pA, qod 


10 


1011011010 


3/8 


0.473404 


pA, qod 


10 


1011011011 


3/8 


0.473404 


pA, qod 


11 


10111111010 


4/9 


0.373716 


pA, qod 


13 


1011011011010 


4/11 


0.479450 


pA, qod 


13 


1011011011011 


4/11 


0.479450 


pA, qod 


15 


101111111111011 


6/13 


0.354176 


pA, qod 


15 


101101111011010 


5/13 


0.467734 


pA, qod 


16 


1011011011011010 


5/14 


0.480804 


pA, qod 



11 



TABLE IV. Periodic orbits extracted from the Be- 
lousov-Zhabotinskii Reaction time-series up to period 15. All 
orbits with a best (normalized) recurrence of less than 0.1 are 
shown. The period 16 orbits are from Ref. [15]. 



Period 



Name 



Recurrence 



Type 



1 



1 



0.016782 



fo 



10 



0.002615 



fo 



101 



0.000128 



fo 



1011 



0.002648 



fo 



10111 
10110 



0.002962 
0.013668 



101110 
101111 



0.006449 
0.029014 



1011110 
1011010 
1011011 



0.005088 
0.041837 
0.010585 



fo 

pA, qod, h 
pA, qod, h 



0.442138 
0.442138 



10111010 
10110111 
10110110 



0.014287 
0.017370 
0.070293 



101111010 
101101011 
101101010 
101101110 
101101111 



0.001720 
0.061843 
0.020884 
0.018380 
0.050094 



pA, qod, h = 0.397081 
pA 
pA 
pA 
pA 



10 
10 
10 
10 
10 

11 
11 
11 
11 

12" 

12 
12 
12 
12 
12 
12 
12 

IT 

13 
13 
13 
13 
13 
13 
13 
13 
14 
14 
14 
14 
14 
15 
15 
15 
15 



1011101010 
1011101011 
1011111010 
1011011010 
1011011011 



0.037910 
0.007237 
0.035966 
0.003440 
0.008971 



fo 
fo 
fo 

pA, qod, h 
pA, qod, h 



0.473404 
0.473404 



10110101010 
10110111010 
10110110111 
10110110110 



0.038044 
0.009119 
0.008447 
0.085593 



fo 



101110101011 
101110101010 
101101010110 
101101111010 
101101101011 
101101101010 
101101101110 
101101101111 



0.031644 
0.013333 
0.029674 
0.055765 
0.032222 
0.016560 
0.048563 
0.064420 



1011110101010 
1011010111010 
1011010101110 
1011011101110 
1011011101010 
1011011111010 
1011011010110 
1011011011010 
1011011011011 



0.052853 
0.022397 
0.004291 
0.000712 
0.011199 
0.046255 
0.097165 
0.013063 
0.004663 



pA 



pA 
pA 

pA, qod, h 
pA, qod, h 



0.479450 
0.479450 



10110101110110 
10110101010110 
10110111010110 
10110110111010 
10110110110110 



0.056515 
0.009366 
0.047037 
0.060024 
0.004920 



fo 



101101011101110 
101101110101010 
101101101011010 
101101101010110 



0.034250 
0.015669 
0.010945 
0.076061 



12 



15 


101101101110110 


0.049261 


pA 


15 


101101101111010 


0.007018 




15 


101101101101010 


0.049573 




15 


101101101101110 


0.010284 




15 


101101101101111 


0.078043 




16 


1011011011101010 


na 


pA 


16 


1011011011111010 


na 


pA 


16 


1011011011011010 


na 


pA, qod, h = 0.480804 


16 


1011011011011011 


na 


pA, qod, h = 0.480804 



TABLE V. Number of periodic orbits extracted and pre- 
dicted for time series data from the Belousov-Zhabotinskii 
reaction (« 500 points in the return map). As illustrated 
with the Rossler data, we expect that the time series is far 
to short to be able to extract all the predicted orbits of any 
except the lowest periods. 



Period 


# predicted 


# found 


# missing 


# marginal 


1 


1 


1 






2 


1 


1 






3 


1 


1 






4 


1 


1 






5 


2 


2 






6 


2 


2 






7 


4 


3 


1 




8 


5 


2 


2 


1 


9 


8 


4 


2 


2 


10 


11 


5 


6 





13 



